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Abstract 

Fermi-LAT has confirmed the excess in cosmic positron fraction observed by PAMELA, which 
could be explained by dark matter annihilating or decaying in the center of the galaxy. Most 
existing models postulate that the dark matter annihilates or decays into final states with two 
or four leptons, which would produce diffuse gamma ray emissions that are in tension with data 
measured by Fermi-LAT. We point out that the tension could be alleviated if the dark matter decays 
into three-body final states with a pair of leptons and a missing particle. Using the goldstino decay 
in a certain class of supersymmetric theories as a prime example, we demonstrate that simultaneous 
fits to the total e"*" + e~ and the fractional e"'"/e~ fluxes from Fermi-LAT and PAMELA could be 
achieved for a 2 TeV parent particle and a 1 TeV missing particle, without being constrained by 
gamma-ray measurements. By studying different effective operators giving rise to the dark matter 
decay, we show that this feature is generic for three-body decaying dark matter containing a missing 
particle. Constraints on the hadronic decay widths from the cosmic anti-proton spectra are also 
discussed. 



I. INTRODUCTION 



The intriguing positron excess in PAMELA [T] has been recently confirmed by Fermi- 
LAT [2] with the help of the Earth's magnetic field to distinguish positrons from electrons 
due to the lack of the magnetic field on board. Fermi-LAT also extends the positron excess 
range up to a higher energy scale, around 200 GeV. In addition, measurements in the 
total + by Fermi-LAT also exhibit some interesting feature well beyond the 100 GeV 
region [3]. These results imply the existence of new sources of primary positrons. The 
new source could be astrophysical, for example, nearby pulsars The more interesting 
possibility is that they may be a sign of the dark matter (DM), due to either the dark 
matter annihilation in the galaxy halo [5], or the decaying dark matter [6]. In the case of 
the dark matter annihilation, it requires a very large boost factor due to the fact that both 
the resulting positron flux and the relic abundance are determined by the same annihilation 
channel. Some of the possible solutions include the Sommerfeld enhancement [7] [8] [9] and 
Breit-Wigner Enhancement (TU] [IT|- On the other hand, the decaying dark matter has no 
such a limitation, but a mechanism responsible for a long lifetime has to be realized in this 
case. 

The possibility that the excesses in the PAMELA positron fraction and Fermi-LAT e~^+e~ 
spectrum come from the dark matter annihilation or decay and the corresponding constraints 
from the antiproton and gamma ray spectra has been the subject of intensive investigations. 
Most of the studies in the literature focused on the two-body and four-body final states of 
the standard model (SM) particles [6], with the latter comes from an intermediate stage of 
two light "portal" particles [8] [12] [13] [11] [15]. The general lesson from these studies is 
that final states have to be dominantly /i's or r's, since a significant fraction of the electron 
final states would produce a sharp feature in the total + e~ flux which is not seen in 
the Fermi-LAT data. However, these charged final states would produce diffuse gamma 
ray emissions through inverse Compton scatterings that are excluded by the Fermi-LAT 
data [IE] [IT] [H] [I9]. There also exist constraints on decays into final state producing 
hadrons from the anti-proton measurements made by PAMELA and on prompt decays into 
photons from Fermi-LAT [20], HESS [21], and VERITAS [22], especially for annihilating 
dark matter. 

In a previous paper [23j , we proposed a novel scenario of the decaying dark matter where 
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the final states include two standard model particles and a heavy missing particle.^ In a 
certain class of supersymmetric theories where supersymmetry (SUSY) is spontaneously bro- 
ken in multiple sequestered sectors, there will be a goldstino in each SUSY breaking sector. 
Only one linear combination of them is eaten and becomes the longitudinal component of 
the gravitino. The others are dubbed as goldstini in Ref . [28] . If the lightest supersymmetric 
particle (LSP) is the gravitino and the next lightest supersymmetric particle (NLSP) is a 
goldstino, and /2-parity is conserved, the goldstino decays through dimension-8 operators to 
three-body final states containing the gravitino plus a pair of standard model (SM) parti- 
cles. As a result, the goldstino could have a long lifetime and be cosmologically stable. For 
suitable SUSY breaking scales in the hidden sectors [23j, the lifetime can naturally be in 
the range of 10^^ — 10^^ sec, which is necessary for the decaying dark matter interpretation 
of the observed e^/e~ excess. Unlike the gravitino, the goldstino interactions with the SM 
fields are not universal and it is easy to come up with scenarios where the goldstino decays 
dominantly to leptons [23], thereby producing the observed excesses. A distinct feature of 
this scenario is that the decay of the dark matter is three-body, with the gravitino escap- 
ing detection and carrying off part of the energy. The SM particle pair in the final states 
has a smooth and soft injection energy spectrum, and consequently, a good fit to both the 
PAMELA e^/e~ and Fermi-LAT e"^ + e~ data can be achieved with a universal coupling to 
all three generations of leptons, which is a welcoming feature from the model-building point 
of view. This is in contrast with the case of two-body or four-body final states where decays 
into e^e~ pair are disfavored. 

In a more general context, one can imagine other models where the dark matter decays 
into three-body final states containing a missing particle. Given the very different energy 
spectra from the well-studied cases of two-body and four-body decays, the scenario of three- 
body decays obviously deserves a detailed investigation to map out the parameter space 
allowed by various cosmic ray measurements such as the antiproton and gamma ray data. 
This is the main purpose of this work. While we use the goldstini scenario as the prime illus- 
trative example, more general setup can be included by studying various effective operators 
giving rise to the decay of the dark matter. We will see that, for the mass range that could 

^ Other possibilities of three-body DM decays were considered in Ref. [Hj , [5S] , [55] , [37j and by A. Ibarra 
et al in Ref. [6]. 
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fit both the e^/e~ and e"*" + e~ excesses, the injection energy spectra are quite insensitive 
to the type of operators mediating the decay as well as the spin of the dark matter. 

This work is organized as follows. In the next section, we first perform fits to PAMELA 
and Fermi-LAT je~ and Fermi-LAT + data for a wide range of the dark matter 
and missing particle masses, using the goldstino decay as an example. Then, we consider 
different types of higher-dimensional operators that lead to three-body decays and compare 
the spectra of the decay products from different decay operators. We also compare their 
fits to the PAMELA and Fermi-LAT /e~ data. There are small variations but they in 



general give similar results. In section III , we study astrophysical constraints on three-body 
decaying dark matter using the gamma-ray and anti-proton data. Our conclusions are drawn 
in section HVl 



II. THREE-BODY DARK MATTER DECAY FOR THE POSITRON EXCESS 

In the scenario of goldstino dark matter [23], the goldstino decays to the gravitino through 
a three-body process. It was shown in Ref. [23] that the leptons pair produced in the 
decay has a smooth spectrum and can give a good fit to the PAMELA positron excess 
and the Fermi-LAT + e~ spectrum, thereby providing a possible explanation of the 
observed anomalies. In this section we perform a general analysis of the three-body decaying 
dark matter explanation of the positron excesses in PAMELA and Fermi-LAT data by 
including the new Fermi-LAT result of the positron excess up to ~ 200 GeV in our fitting. 
In particular, we vary both the dark matter and missing particle masses to obtain a best fit 
and go beyond the goldstino scenario by studying various effective operators mediating the 
three-body decay. 



A. Goldstini 



Goldstini arise when there are multiple sequestered sectors which break SUSY. For sim- 
plicity, let us consider that SUSY is spontaneously broken in two hidden sectors, then there 
is a goldstino in each sector. If the superpartners of the SM particles receive SUSY breaking 
masses from both sectors, there will be couplings of the SM particles and their superpartners 
to each goldstino. If the SM superpartners are heavier than the goldstini, upon integrating 
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out the superpartners, we will obtain dimension-8 operators between a pair of goldstini and 
a pair of SM particles.^ One linear combination of the goldstini is eaten and becomes the 
longitudinal mode of the gravitino. The other obtains a mass due to supergravity effects 
[28] . If there is a hierarchy in the SUSY breaking scales of the two hidden sectors, the 
eaten goldstino mostly comes from the sector with a larger SUSY breaking scale. The un- 
eaten goldstino, which we assume to be the dark matter, is made mostly of the goldstino 
in the sector with a smaller SUSY breaking scale. After going to the mass eigenstates, the 
dimension-8 operators contains interactions which allow the uneaten goldstino decays to a 
gravitino and a pair of SM particles. If -R-parity is conserved, the lifetime is naturally longer 
than the age of the universe and, with suitable choices of SUSY breaking scales of the hidden 
sectors, can be the required time scale to explain the e+/e^ excess observed by PAMELA 
and Fermi-LAT. 

The goldstino decay operators were derived in the previous paper [23]. We have for SM 
fermions, 

Ml) 1 / m? tan 6^ - mo cot 6' \ „ ,- , , 

4/ = — d.md^iGLq) + h. c. , (1) 

Jeff \ J 

where ( is the goldstion, G^, is the logitudinal mode of the gravitino, q is the SM fermion, 
feff = a/ + /I is the effective total SUSY breaking scale with /i 2 being the SUSY breaking 
F-terms of the two SUSY breaking sectors, tan 6* = /2//1, ?Tii2 are the soft SUSY breaking 
mass of the superpartner of q coming from the two SUSY breaking sectors respectively, and 
m| = ml + m^. For gauge bosons, 

Jeff \ "^A / 

where F{F) = F^j^^a^^ {F^ylf^^) is the gauge field strength tensor, mi_2 are the SUSY breaking 
gaugino masses coming from the two SUSY breaking sectors, and mx = rhi + m2. For the 

^ The operators involving the Higgs fields actually have a lower dimension. However, it is suppressed by 
the same SUSY breaking scale with the mass dimension made up by the Higgs mass parameters in the 
numerator. 
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Higgs field, there are terms from botli Kaliler potential and superpotential, 



(1) ^ 

2H ^,f2 



+u ^ d + h.c. , (3) 



eS 



+u ^ d + h.c. , (4) 



where 



k = > = mi„tan6' - m2„cot6' , a = u,d, (5) 

i 

B = , 5B = Bi tan - B2 cot ^ , (6) 

j 

with i = 1, 2 indicating the SUSY breaking sector where the soft SUSY breaking parameters 
come from. The operators with the Higgs fields can actually induce two-body decay ( — )■ 
Gih after substituting the Higgs vacuum expectation values. The two-body decay mode is 
suppressed by a factor (f /m^)^. Compared with the phase space suppression of the three- 
body decay, the two-body decay rate is expected to be of the same order as the three-body 
decay rate but may be somewhat larger depending on the model parameters. 

Unlike the gravitino, couplings of the goldstino to SM fields are not universal, but rather 
depend on the fractions of the soft SUSY breaking parameters of the corresponding super- 
partners coming from the two SUSY breaking sectors. To explain the cosmic e~^/e~ excess 
while satisfying other cosmic ray constraints, the decays should dominantly go to leptons. 



The constraints on the branching fractions to other particles will be studied in Sec. [HT 

To the leading order, the goldstino acquires a mass which is twice that of gravitino due 
to supergravity effects [28]. However, the relation can be modified significantly beyond 
the leading order [29]. In this subsection, we extend the analysis to explore a wider range 
of the goldstino mass and the gravitino masse '^g^ by relaxing the mass relation of 
= 2mg^. We follow exactly the same procedure as in Ref. [23j. To recap, we make 
use of the Bessel function method of Ref. [30J to obtain the positron flux detected on the 
Earth due to the DM decay by assuming the MED model parameters. For the background 
fluxes we employ the "model" presented by the Fermi-LAT collaboration [31], which can be 
analytically parameterized as shown in Ref. [32]. The Moore profile in Ref. ^3] is used for 
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FIG. 1: The best fits to the e~^/{e + e"*") ratios observed by PAMELA and Fermi-LAT, and the 
Fermi-LAT + e"*" spectrum for several goldstino masses. 

the DM halo. With the assumption that the dark matter density is = 0.3 GeV/cm^, we 
perform combined fits to both PAMELA and Fermi-LAT data, including the new Fermi-LAT 
positron fraction data [2] , by varying the decaying DM lifetime and the overall normalization 
of the primary component of the background flux, as described in Ref. [32] • Because the 
e^/e^ flux with energies below 10 GeV measured at the top of the atmosphere is significantly 
influenced by the solar modulation effect, we only include data points above 10 GeV from 
PAMELA for the total x^- We assume universal couplings to all three lepton flavors which 
give a better fit than couplings to any single flavor. 

The fits to the PAMELA and Fermi-LAT data for several sample goldstino masses are 
shown in Fig. [ij The P^r degree of freedom (d.o.f ) of the fits as functions of the goldstino 
mass m,^ and the gravitino mass are shown in Fig. [2j The solid line in Fig. 



2 a 



is the 



best fit x^/d.o.f. for a given m,^ by varying ttl^^ freely, only subject to ttl^^ < m^. We see 
that x^/d.o.f. asymptotically approaches < 1 when > 2 TeV. The dashed line, however, 



assumes = 2mg^, in which case the minimum of occurs around m^^ ~ 2 TeV. Fig. 2(b) 
shows the contour plot of versus — mg^ , from which we see that ~ 1 occurrs when 

— mg^ ~ 1 TeV and > 2 TeV. Consequently, we will use = 2 TeV as a benchmark 
value for later discussion. 

The underlying reason for these features is the hardening feature around 300 ~ 500 
GeV in the e"*" + e~ flux from Fermi-LAT, which dictates — ttl^^ to be around 1 TeV. 
This number comes about because we assume universal couplings to all three lepton flavors 
and the resulting softer energy spectra due to the fj, and r components drives the value of 




oou' ■ ■ ■ — ■ — ■ ■ ■ — ■ — ■ — ' ' ' ' 

10 10 10 10 10 

(GeV) ,77.^ (GeV) 



(a) (b) 

FIG. 2: These plots display relations between x^/d.o.f. of best fits to both PAMELA and Fermi- 
LAT data and the goldstino mass m^. On the left plot, we vary and fn^^ independently (solid 
line) as well as dependently by assuming the mass relation: = 2mg^. The right plot shows 
contours on the and rrif^ — inn^^ plane. 

m(^ — m^^ to be larger than twice the energy scale at which the hardening spectrum appears: 
— rn^^ > 600 — 1000 GeV. On the other hand, the goldstino mass m,^, which sets the 
overall energy scale, need to be heavier than 2 TeV since a smaller mass will produce a 
feature in the total e"^ + e~ flux at around 500 — 1000 GeV, which worsens the fit with the 
Fermi-LAT data. (See Fig. |l|) 




mg^ (GeV) 

FIG. 3: The contour plot represents 90%, 95% and 99% C.L. for fits to PAMELA and Fermi-LAT 
data by varying the rrip and the lifetime r of C,. 
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In addition to varying masses of the goldstino and gravitino, we also consider the best 
fit hfetime r of the decaying goldstino, which is determined by the strength of the excess 
in the PAMELA /e~ flux in the energy regime of 10 — 200 GeV. We present a contour 
plot corresponding to 90%, 95% and 99% C.L. of fits to PAMELA and Fermi-LAT data in 
the versus r plane in Fig. |3j It shows that for mc^ = 2 TeV, the best fit corresponds 
to fn^^ ~ 900 GeV. As fn^^ increases (decreases), it results in the softer (harder) injection 
spectrum, which can be compensated by a shorter (longer) lifetime. 

Before switching to three-body kinematics, we would like to comment that for > 2 
TeV, the superpartners of the SM particles need to be heavier than a few TeV which implies 
some fine-tuning of the electroweak breaking scale. However, the absence of the SUSY signal 
at the LHC so far indicates that the superpartners may indeed be heavier than expected if 
SUSY exists at the TeV scale. The hint of a Higgs boson around 125 GeV is also consistent 
with a heavy SUSY spectrum. 

B. Three-Body Dark Matter Decay Kinematics 

If a dark matter particle decays to two-body final states, the energies of the decay prod- 
ucts in the rest frame of the DM particle are fixed. In contrast, the three-body decay will 
produce softer and broader spectra as we saw in the goldstino example. One can imagine 
there are other models where DM particles have similar three-body decays into a missing 
particle and two SM particles. The spectra of the decay products depend on the properties 
of DM particles and their interactions, but do share some generic features. By parametrizing 
different three-body decay mechanisms by effective higher dimensional operators, we will ex- 
amine the similarities and differences of various three-body decay spectra in this subsection. 

We start by comparing the injection energy spectrum of the SM fermions from decays via 
dimension-six four-fermi interactions with that from the benchmark scenario of goldstino 
decays, which occur through the dimension-eight operator in Eq. ([T]). More specifically, 
we consider the following Dirac structures |27j: (pseudo-)scalar, (pseudo-)vector and tensor 
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interactions, which are written as 

Ai 

A 
A 

A 
A3 

A4 

A5 



^2^X*DM^g*g + h.C., 
-^Wx7^^DA/^g7^*<7 + h. c, 

^2^x7^^DA/'*g7/.^g + h. c, (7) 
^2^X7''7^^DAf^g7/.7^^<? + h. c, 



— ^XfT^'^^DA/^gCT^i.^g + h. C, 

where we assume both the dark matter (DM) and the missing particle in the decay product 
(X) are Dirac fermions, while q refers to the SM fermions. The decay process induced by 
these operators is 

DM ^ X + g + g . (8) 



Fig. 4(a) shows the comparison of energy spectra of SM fermions from the three-body 
decaying dark matter, where the total width is normalized to unity and the masses are 
fixed to mDM = 2mx = 2 TeV. The goldstino decay results in the hardest spectrum due 
to fact that goldstino is derivatively coupled, which gives more weight to the higher energy 
modes. However, one sees that all spectra exhibit similar qualitative features, peaking in 
the 400 — 600 GeV region. Therefore, we expect all of them to give rise to similar fits to the 
Fermi-LAT and PAMELA data points. 

In addition to decays to SM fermions, we also compare various decay operators into SM 
gauge bosons, as they will be constrained by the anti-proton data from PAMELA as well as 
the gamma ray data from Fermi-LAT. Here we list three types of different operators to be 
compared with the benchmark scenario of goldstino decay in Eq. (|2]): 



Ah/f^ 



^3 ^x'^DMW;W''^ + \i.C., 

^-^x^DMF^^F'^^^ + h.c, (9) 

^Wx7^v^*dmf;.f;. + h. c, 

where F^j, can be any of the SU{?>), SU{2), and U{1) SM gauge bosons. The suppres- 
sion of f ^ in front of the first operator signals its origin from the gauge invariant operator 



{\h' / ^?)^ DA[D^ D fj_H . We see in Fig. 4(b) that, again, various operators have very 



similar qualitative features, even more so than the energy spectra of the SM fermions. 

10 



DM=2 Tc-V X=l Tc-V 



DM=2 Tc-V X=l Tc-V 













Gloldstino dcciiy 










if' / 








#x7''7^*Dfl/*q7^7'^*y/A^ 






#x7^*'DM*,7^*g/A^ 






f^o-'"'4'£)M4'giT^^*^/A^ 





— Goldstino decay 



100 200 300 400 500 600 

Energy (GeV) 




100 200 300 400 500 600 700 800 

Energy (GeV) 



(a) (b) 

FIG. 4: The comparison among different DM decay mechanisms into fermions and WW, by nor- 
malizing the total decay width F = 1. We assume the mass of the dark matter is 2 TeV and 1 TeV 
for X and Gl- 
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FIG. 5: The comparison between spin-0 and spin-l/2 dark matter with several decay mechanisms. 
We assume moM = 2mx = 2 TeV. 



For the higher dimensional operators involving a pair of Higgs fields, if the Higgs fields 
are not both derivatively coupled, substituting in the Higgs VEV will induce a two-body 
decay DM hX. The Higgs boson will have a delta-function spectrum independent of 
the exact operators. We shall see that this two-body decay gives comparable constraints to 



those from the corresponding three-body decay in section HI 



We conclude this subsection by comparing energy spectra of three-body decaying scalar 



dark matter with those of fermionic dark matter. From Fig. 5(a) it is notable that, for 
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FIG. 6: Fits to PAMELA and Fermi-LAT using two-body decaying dark matter. For the upper 
panel we include data on both positron fraction and total e~ + e"*" flux in the fit, while in the lower 
panel we only include the positron fraction data in computing the ■ 



mx = m£)j\//2, spectra of scalar dark matter follows closely the counterparts of fermionic 
dark matter. However this behavior disappears in the limit of massless missing particle mass, 
nix = 0, as shown in Fig. 5(b)[ Therefore, in terms of fits to PAMELA and Fermi-LAT, we 
expect to have similar x^/d.o.f. for both scalar and fermionic dark matter. 



C. e+ and e 



Having established that the injection energy spectra of three-body decaying dark matter 
do not depend sensitively on the dynamics giving rise to the decay in the mass range of 
interest, we now illustrate how three-body decaying dark matter could give better fits to both 
PAMELA and Fermi-LAT /e~ data than the conventional model of two-bdoy decaying 
dark matter. 



12 



2 TcV Goldstino 
900.S GeV Gravitino 
r = 1.44lS:»»xlO" 




— Goldstino decay + Background 

— Background 
PAMELA 
FERMI 



10 10 

Energy (GeV) 

DM ^ Xl+l- (*A'7''*m/*,7/I'y 



2 TcV DM 
(i34.1 GcV X 




DM decay + Background 

Background 

• PAMELA 

• FERMI 

10^ 



Energy (GeV) 
DM ^ Xl+l- (*x7'*DM*,7''*,: 



2 TeV DM 
766.0 GeV X 




- DM decay + Background 

— Background 
PAMELA 
FERMI 



10 10 

Energy (GeV) 



2 ToV Goldstino 
900.8 GeV Gravitino 
r = 1 44+" ''^ xlO^^ sc 



- CSoldstino decay + Background 

- Background 
FERMI 




10 10 10 

Energy (GeV) 

DM ^ Xl+l- (f x7^*DA/*,7^*,) 



2 TeV DM 

6;U.l GeVX 

T = 1.3^)^^l^'" xlO^^ sec 



- DM decay + Background 

- Background 
FERMI 




10 

Energy (GeV) 
DM Xl+l- (*x7^*DM4',7'^*q) 



2 TeV DM 
7fjfi.O GeV X 



DM decay + Background 

Background 

• FERl\n 




10 

Energy (GeV) 



FIG. 7: The fits to PAMELA and Fermi-LAT between different operators. Three-body decays 
have the similar o-nd lifetime, controls the hardness of injection energy spectra; therefore, 

goldstino has the largest to balance the hardest spectrum. 



We first examine the fits for two-body decays of dark matter to a pair of leptons witli 
universal couplings to all three generations. As explained before, the observed excess of 
positron fraction from both PAMELA and Fermi-LAT can be accounted for by varying the 
lifetime of the decaying dark matter as long as the injection energy spectrum is not too 
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soft. However, in order to explain the hardening feature in the total e~ + flux observed 
by Fermi-LAT, the injection spectrum must peak around (9(400) GeV. As a consequence, 
the two-body decay has difficulty fitting both the positron fraction and total flux, since 
the injection energy is a delta-function peak at half the mass of the dark matter in the 
centre-of-mass frame, which is demonstrated in the upper panels of Fig. [6j Note that one 
could use only the PAMELA and Fermi-LAT /{e~ +e^) data in computing the x^/d.o.f., 
then the positron excess can be well described with a shorter lifetime compared to those of 
three-body decay cases. But the resulting total (e~ + e"*") flux is in severe tension with the 
Fermi-LAT data, as can be seen in the lower panel of Fig. [6] 

The best fits for several three-body decay mechanisms are displayed in Fig. [7| again 
assuming universal couplings to all three lepton flavors. We see the overall fits to both 
PAMELA and the Fermi-LAT are better than the case of two-body decay. Among the 
different decay mechanisms, the goldstino case has largest missing particle mass, which 
softens the harder energy spectrum resulting from its derivative couplings. Similarly, the 
pseudo-scalar operator needs a larger mx than the vector one due to a harder injection 
spectrum, but a shorter lifetime, driven mostly by the excess of e^/ (e^ +e^), to compensate 



for the flatness in the spectrum shown in Fig. 4 (a' 



So far we have assumed universal couplings to all three lepton flavors. However, the 
(pseu do-) scalar operators break the chiral symmetry and one may expect that the coefficients 
are proportional to the fermion masses, which implies that the decay to r leptons will be 
dominant. In this case, these two operators can not explain both PAMELA and Fermi-LAT 
data. The spectrum of e+(e~) from r decays is quite soft and would populate mostly the 
low-energy region. In order to explain the excess in PAMELA data, the lifetime has to be 
increased significantly and the missing particle mass mx has to be reduced, which results in 
a much harder spectrum than the Fermi-LAT total e~ -|- e"*" flux data as shown in Fig. ^ 



III. ASTROPHYSICAL CONSTRAINTS 



In this section we consider astrophysical constraints on three-body decaying dark mat- 
ter. There are three main categories: 1) the diffuse gamma-ray due to inverse Compton 
scatterings (ICS) and final state radiation (FSR), 2) prompt photons that are direct decay 
product of the dark matter, and 3) anti-proton flux measurements which constrain decays 
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FIG. 8: Best fits for the pseudo-scalar operator to PAMELA and Fermi-LAT combined, assuming 
couplings proportional to the fermion mass. The PAMELA (Fermi-LAT) excess can be accounted 
for properly at the cost of a longer lifetime and mx = 0, which yields an unwanted bump in the 
Fermi-LAT e~ + e'^ data. 

into hadronic final states. 



A. Diffuse 7 

It is known that the leptonic final states from dark matter decays yield photons via FSR 
as well as ICS when the produced leptons interact with background photons. Both are 
subject to the constraints from the Fermi-LAT gamma ray data. In this subsection, we 
study such constraints by assuming the e^/e~ excess in PAMELA and Fermi-LAT is the 
consequence of dark matter decays. Since we have shown that different three-body decay 
operators have similar best-fits to e+/e^ data for 2 TeV dark matter mass, it is sufficient to 
focus on the goldstino case only. 

We make use of the Fermi-LAT gamma ray data in Ref . [M] , which provides information 
of diffuse Galactic emission (DGE) and extragalactic gamma ray background (EGB). For 
DGE we fit to the "Galactic diffuse (fit)" data from Table I in Ref. [34j, which is DGE 
averaged over the Galactic latitude range \b\ > 10° as measured by Fermi-LAT. The DM 
signal results from the sum of ICS and FSR produced by charged leptons, the products of 
the goldstino decay. The averaged photon fiux from FSR over a solid angle of interest from 
Table 3 in Ref. [35] and a numerical code for ICS in Ref. [36j are utilized for calculating the 
photon fiux from the Galactic dark matter decay. For the fitting procedure, we vary the 
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normalization of the background, whose shape is taken from the "Galactic diffuse (model)" 
numbers in Table 1 of Ref . [31] , to minimize the for various goldstino decay widths. In 
addition, for the goldstino decay, the gravitino mass is extracted from the best fits to the 
/e~ and + data. 

Before presenting the results, we would like to comment on the effect of different electron 
propagation models, dubbed MIN, MED, and MAX [30] [38]. For FSR, secondary pho- 
tons are produced near where primary leptons are produced and once produced, photons 
propagate directly toward the Earth; therefore, they are independent of the propagation 
models of electrons through the galaxy. On the other hand, ICS photons are produced from 
electrons scattering off background photons while the electrons propagate and lose energy. 
Hence, it has dependence on the propagation models of electrons. For the MIN model, the 
low/intermediate energy electron (< 200 GeV) will be more easily stopped and fail to pro- 
duce ICS photons. Fermi-LAT gamma ray data ranging from 0.2 to 100 GeV are influenced 
by this effect. Therefore, for both two- and three-body decay, decay scenarios adopting the 
MIN propagation model are less constrained compared to MAX and MED. MAX has an 
opposite feature, i.e., producing more low/intermediate ICS photons but the effect is less 
dramatic than MIN. In other words, MAX is similar to MED. 

We present results in Fig. [9] corresponding to the 90% and 95% C.L. upper (lower) limit 
on the decay width (or equivalently lifetime) of the goldstino, normalized to T^^'^^ from best 
fits to e"*" and e~ data, based on the Fermi-LAT DGE data. Again, the photon flux resulted 
from the goldstino decay consists of both ICS and FSR contributions. The goldstino decay 
can satisfy the diffuse gamma ray constraints. In contrast, for two-body decays with the 
MAX or MED propagation model, the decaying dark matter is ruled out as an explanation 
to the excess of and e~. In the case of MIN, both two- and three-body decay can avoid 
the gamma ray constraint as explained before. Note that FSR contribution is numerically 
subdominant compared to those of ICS on both two- and three-body cases. The reason is 
that the r decay channel, which produces relatively more photons via showering than e and 
/i channel, gets diluted by the assumption of universal leptonic couplings. Therefore, if the 
dark matter decays into r only, FSR would be far more stringent than ICS as demonstrated 
in Refs. [HI [18]. 

On the other hand, the Fermi-LAT EBG data provide constraints on the isotropic diffuse 
gamma ray. From Ref. [IS], the DM decay contributes to the isotropic flux in two ways. 
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FIG. 9: The goldstino decay can satisfy the Fermi-LAT DGE gamma ray constraints. In contrast, 

for the two-body decay, only MIN is allowed by virtue of the reason mentioned in the text. Note that 

for the certain range of the DM lifetime, the existence of the DM signal compensates the difference 

between the data and background and in turn results in dips in fits. 

which can be expressed in terms of the differential flux, 

'j^'^ Isotropic 



+ 4tt 



dE^dQ 



(10) 

where the first term is an isotropic extragalactic cosmological flux from the decays at all past 
redshifts and the second term is the residual contribution from the Galactic DM halo. Note 
that the latter is of course not isotropic but the minimum will be the irreducible component 
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of the isotropic gamma ray flux. We employ the same assumption as in Ref. [19] that the 
minimum of Galactic contribution is located at the anti-Galactic Center. With the help of 
the code, EGgammaFluxDec in Ref. which is based on the two-body decay and includes 
both FSR and ICS from the primary charge leptons ^37], we properly convolute it with 
our three-body injection spectra to obtain the isotropic cosmological flux. For the fltting 
procedure, we again use the EGB data from Table I of Ref. [34], which can be described 
very well by a featureless power law with index 7 = 2.41 ± 0.05, and hence a power law 
background. The normalization of the background and the index of the power law are being 
varied to flnd the minimum of for different goldstino lifetimes. The results are presented 



in Fig. [10| Based on the numerical results, the cosmological flux, which is independent of the 
propagation model, is dominant over that of the DGE minimum and clearly the bounds on 
the goldstino lifetime are similar for different propagation models. Therefore, the EGB data 
impose more stringent constraints than those of DGE, i.e.. Fig. |9j subject to uncertainties 
on the propagation models. We can see that all of two-body decay are excluded and even 
the goldstino decay in the case of MIN is also disfavored, since the lifetime (decay width) of 
the goldstino is shorter (larger) in MIN to account for fewer electrons and positrons reaching 
the Earth. 



B. Prompt 7 

In this subsection, we explore the situation where photons are direct products of three- 
body dark matter decays, i.e., primary photons.'^ In this situation, we have two independent 
parameters: the decay widths of goldstino into leptons and into photons. Moreover, the 
results are obviously independent of the electron propagation models like FSR and the 
cosmological diffuse flux, and the constraints would be model-dependent. As the injection 



spectra of different operators from II B are similar, we expect those four operators to have 



the similar results. Here, we also involve both DGE and EBG data with the same fltting 



procedure described before. From Fig. 11, the constraint mainly comes from EGB data 
and the goldstino decay can not be ruled out even if the goldstino decays 80% into photons 
and 20% into leptons, which accounts for the e^/e~ excess observed by PAMELA and 



^ Note that the signal of the two-body prompt decay is situated beyond the reach of the Fermi-LAT gamma 
ray data for 2 TeV dark matter |34) . 
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FIG. 10: The constraints from the Fermi-LAT EGB data. Only the goldstino decay with MAX and 
MED pass the test of the isotropic diffuse gamma ray flux. Like in Fig. we see the dip appears 
especially for the goldstino decay, but here the power-law index of the background is being varied 
and therefore the shape of the background changes for the different DM lifetime. In addition, the 
change on the power-law index is also the underlying reason for the local minimum around Ax^ ~ 5 
in the case of the two-body decay in MAX and MED. 

Fermi-LAT. The underlying reason is, for 2 TeV goldstino, the prompt photons are mainly 
located beyond 100 GeV and consequently the DGE data yield a very loose bound. In these 
cases. Air Cherenkov Telescopes such as Veritas are well poised for DM decay to photons 
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FIG. 11: The Fermi-LAT DGE constraints (upper panels) and EGB constraints (lower panels). 
The 2 TeV goldstino decay to photons could easily pass the test of the Fermi-LAT gamma ray data 
because the injection spectrum of photons peaks well beyond 100 GeV. On the other hand, a 500 
GeV goldstino is severely constrained by Fermi gamma ray data. 



of O(TeV) energy. On the other hand, the cosmological diffuse photons will get redshifted 
and constrained by the EGB data. For comparison, 500 GeV goldstino would be severely 



constrained by both of DGE and EGB data as shown in Fig. 11 



C. p 

In addition to measurements of the spectra made by PAMELA and Fermi-LAT, the p 
spectra has also been well measured by PAMELA. While the species show an excess at 
high energies, the p measurement shows a spectra consistent with astrophysical sources |39] . 
The most recent PAMELA anti-proton measurement from 60 MeV to 180 GeV has further 
confirmed this result jlO]. This agreement with the expected background can serve to limit 
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FIG. 12: Antiproton and electron energy distribution functions of fi,T-leptons, u-quarks, gluons, 
W, Z and Higgs boson decay. The hardness of the electron in the lepton and W/Z decays result 
from prompt production, whereas the quark/ gluon decays are very soft, producing electrons only after 
QCD showering and hadronization. Similarly, the antiproton energies are very soft and are similar 
for quarks and W bosons, with a slightly harder component for W decays. For two body decay, 
the resulting positron and antiproton spectra are similar to these distributions. The corresponding 
three body decay convolves the intermediate state energy with these distributions. 



the total annihilation rate of dark matter to hadronic final states [H] as well as hadronic 
interactions of the dark matter, which could have important implications for direct detection 
experiments [12]. 

We study the restrictions of hadronic decay modes of the dark matter by the potential 
contribution from the decay to qq, gg, W~^W~ , ZZ and hh modes, with nih = 125 GeV. In 



the left panel of Fig. [12] we show the fragmentation function of the above final states into 
an anti-proton. Since these final states may also produce electrons upon further decays, it 
is important to study whether the partial width allowed by the anti-proton spectra could 
result in additional noticeable contributions in the positron measurements that were not 



included in our fits in the previous section. Thus in the right panel of Fig. 12 we also 
show the fragmentation function of all possible final states into the positron. If the dark 



matter decays through two-body kinematics, then Fig. 12 gives precisely the injection energy 
spectra of the anti-proton and the positron. For three body decay, these energy distributions 



21 



are convolved with the intermediate state's energy distribution. It is clear that due to the 
nature of fragmentation and hadronization, all non-leptonic modes provide a soft anti-proton 
spectra. 

Once produced, the anti-protons propagate through the galaxy under the effects of dif- 
fusion due to the galactic magnetic field, the convective wind away from the plane of the 
galaxy and annihilations with interstellar protons; we model the anti-proton propagation 
according to Ref. |13]. As with the modeling of the e"*" propagation, we primarily adopt 
the MED model and assume the dark matter is distributed according to the Moore profile. 
In addition, once the anti-protons approach Earth, solar modulation effects alter their low 
energy spectra. We include this effect for anti-protons since the energy range for solar mod- 
ulation is well positioned within the PAMELA data. To this end, we adopt a Fisk potential 
of = 500 MV [B]. 

To illustrate how consistent the two-body and three-body decay scenarios are with the 
anti-proton data, we fit the prediction of each model to the PAMELA p/p data [10]. We 



take the best-fits to the electron data in Sect. II C and simultaneously vary the dark matter 
signal normalization and total p background normalization. We show in detail the fits for 
the qq and W^W~ mode in the MED propagation model. We later summarize the results 
for other decay modes and propagation models. 

The two-body decay of dark matter yields two particles with well defined energy. If they 



are massless, the resulting anti-proton spectra appears identical to Fig. 12 with Eparent = 
mnM/'^- The energy of injected anti-protons from massive parents are altered with two- 
body kinematics. In Fig. [13} we show the two-body dark matter decay to quark modes 
(left panels) and W boson modes (right panels) for moM = 2 TeV. The comparison with 
the PAMELA anti-proton data (top panels) shows that a considerable anti-proton flux is 
possible in energies beyond the present data. The Ax^ fits (middle panels) show no preferred 
value for the partial width, as expected since no excess is observed in this data. Had there 
been a preferred decay rate, the position of the minimum would not be asymptotically 
approaching zero. Finally, we show in the bottom panels the components of the proton 
and anti-proton flux based on the 95% CL. flt. Overall, we see that decay rate to the 
quark mode at 95% CL. must be smaller than ^ 1/(2 A x 10^^ s) and the rate for the 
W^-boson mode must be smaller than ^ 1/(1-3 x 10^'' s). One could also compare these 
bounds on the hadronic decay width with the best flt decay widths into the lepton modes 
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from fitting tlie electron/positron data. If we take the best-fit lifetime from fitting two- 
body decays to both the positron fraction and the total e+ + e~ flux, which is shown in the 
upper panels in Fig. [6| the hadronic partial widths fall around C(10%) of the lepton modes. 
The upper limit on the two body decay of dark matter in the VT-boson mode is roughly 
Tdm^w+w- < 0.17 X Tdm^£+£- from the PAMELA anti-proton data. Likewise the limit 
on the q mode is TuM^qq < 0.09 x TuM^e+e-- However, these ratios are a factor of 2 — 3 
smaller if one only includes the positron fraction in the fit for the leptonic widths, which 
was demonstrated in the lower panels of Fig. |6j 

The quark mode is more constrained by roughly a factor of two relative to the WW mode 
since the anti-protons originate from quark hadronization. In comparison, the ly-boson can 
decay to iu, giving a reduction in the overall p rate. Moreover, as the VT-boson decays 
to pairs of quarks, the resulting anti-protons are softer compared with the dark matter 
decay to quarks. However, below Tp = mpEi^j/mw, the p spectra is cut off. This can be 
understood since the p's are predominantly relativistic, meaning they have a kinetic energy 
of at least 0{mp) in the W boson rest frame. Boosting to the dark matter rest frame places 
the energy cutoff near Tp = mpE\^^/mw- Overall, the flux from lower energy anti-protons 
becomes buried in the increasing background p flux, thus providing a lower p/p contribution. 
This is clearly seen in the position of the peak of the DM signal in the lower panels. The 
propagation model dependence is shown in Table |Tj The decay lifetime fit to the lepton 
modes from the PAMELA data are shown in units of 10^^ s. The subsequent 95% CL. 
lower limits from the PAMELA p data on the decay lifetime of the two body decay modes to 
gg, gg, W~^W~ , ZZ, and hh are also summarized. In parentheses are the respective ratios 
between the leptonic lifetime fit and the hadronic lifetime limit. We see that obviously the 
most constrained modes are the purely hadronic qq and gg modes. Generally, within the 
MAX propagation model, more suppression into the decay of hadronic final states is required 
to maintain agreement with data. Within the MIN propagation model, the hadronic states 
have to be suppressed no more than an order of magnitude. 

Generically, the three-body decays through various operators of dark matter to quarks 
and ly-bosons create similar p spectra. We therefore concentrate on the analysis of the p flux 
in the goldstino decay scenario where it decays through the dimension-8 operators. Since 
the injected q and W in the three-body decay is of lower energy, the resulting anti-protons 
are softer than the two-body case. As a result, the constraint on the allowed decay width in 
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FIG. 13: Contribution to the p/p ratio from dark matter 2-hody decay for the W and q decay 
modes. Since the W -boson tends to produce fewer p than the q-mode, the decay rate allowed by the 
PAMELA p data is larger by nearly a factor of two. Dark matter decay rates for W-bosons (q) as 
high as 0.17(0.09) x r£)^,/_^£+^- can be accommodated by the PAMELA p data. 
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FIG. 14: Similar to Fig. \13[ but for dark matter three-body decay via the Dim-8 goldstino decay 
operators. Decay rate for W-bosons (q) as high as 0.14 (0.11) x ^ Df^.f^Q^i+i- can be accommodated 
by the PAMELA p data. 
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TABLE I: Dark matter two-body decay mode lifetime dependence on the propagation model for fits 
to the PAMELA/Fermi e+, Fermi e+/(e+ + e~), and PAMELA p data with = 2 TeV. The 
decay to the i^i^ mode is the best lifetime fit provided by the Fermi-LAT and PAMELA electron 
data, while the other decay modes are 95% C.L. lower limits on the lifetime and the ratio with 
respect to the best fit lepton lifetime in parenthesis. 



lO^e s 


C ^ i+£- 


C^qq C^gg C ^ w+w- zz C^hh 


MAX 

MED 
MIN 


2.17 
2.13 
2.01 


49.7 (22.9) 80.9 (37.3) 29.6 (13.7) 31.4 (14.5) 47.5 (21.9) 

23.8 (11.2) 39.2 (18.4) 12.8 (6.01) 13.7 (6.42) 21.0 (9.86) 
5.63 (2.80) 9.30 (4.62) 2.89 (1.44) 3.09 (1.54) 4.78 (2.38) 



TABLE II: Goldstino decay mode lifetime dependence on the propagation model for fits to the 
PAMELA/Fermi e+, Fermi e+/(e+ + e~), and PAMELA p data with = 2 TeV. The decay to 
the GL£^i~ mode is the best lifetime fit provided by the Fermi and PAMELA electron data, while 
the other decay modes are 95% C.L. lower limits on the lifetime and the ratio with respect to the 
best fit lepton lifetime in parenthesis. 



10'^>^ s 


C ^ GLi+r 


C ^ GLqq C ^ GLgg C ^ GlW+W- C ^ GlZZ ( ^ G^/i/i C ^ G^/i 


MAX 
MED 
MIN 


1.51 
1.44 
0.89 


23.2 (15.3) 23.6 (15.6) 22.8 (15.1) 21.8 (14.4) 22.0 (14.5) 23.8 (15.7) 
12.5 (8.68) 14.7 (10.2) 10.6 (7.36) 10.2 (7.08) 10.0 (6.94) 10.5 (7.29) 
2.48 (2.83) 2.89 (3.29) 2.10 (2.39) 1.99 (2.27) 1.99 (2.27) 2.39 (2.73) 



every hadronic channel is weaker than the corresponding two-body case, which can be seen 
by comparing Table [ll] with Table [l} The decay rate to the quark mode at 95% C.L. must 
be smaller than ^ 1/(1-3 x 10^^ s) while the rate for the VT-boson mode must be smaller 
than ^ 1/(1-1 X 10^"^ s). In terms of the ratio of the hadronic decay rates over the best-fit 
leptonic decay rates, the 95% C.L. upper limit on the three-body decay of dark matter in the 
VT-boson mode is roughly Tdm^w+w- < 0.14 x T^M-^e+i- from the PAMELA anti-proton 
data. Likewise the limit on the q mode is TriM^qq < 0.11 x Tdm^i+i-- Similar to the Table[T] 
for the two-body decays, the propagation model dependence for three-body decays is shown 
in Table M 
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FIG. 15: The effect on fits to PAMELA and Fermi-LAT from inclusion of hardonic final states. 
The bounds is much looser than ones from antiproton data which justify our fitting procedure. 



It is worthwhile to mention that, for the three-body decay, in principle, we should simul- 
taneously fit both of the positron (electron) and the antiproton data by including all kinds 
of decay products such as leptons, quarks and gauge bosons, etc. Instead, we fit positron 
data first with universal couplings to leptons only and obtain or mx, and then fit 

antiproton data with the extracted mg^(or nix) to see how much the hadronic final state 
is allowed. In this procedure, first it is easier to pin down the sufficient conditions for good 
fits to positron data and second it is numerically faster on searching for minima. Neverthe- 
less, one can argue that including other final states might change best fits to positron data 
significantly since these final states yield positrons (electrons) as well, and therefore modify 
the whole picture. Fig. [15] clearly shows that inclusion of hadron final states will have a 
modest effect in terms of on fits to PAMELA and Fermi-LAT positron data even if the 
quarks or VT-bosons comprise more than 25% of decay products, which is much higher than 
the upper bounds from the antiproton data consideration. Similar conclusions hold for the 
effect on the gravitino mass and the goldstino lifetime. As a result, our way of fitting is 
properly justified. 



IV. CONCLUSIONS 



In this work we have attempted to make a strong case for the three-body decaying dark 
matter, by studying how the three-body kinematics could allow simultaneous fits to both the 
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PAMELA and Fermi-LAT positron excesses without being excluded by other astrophysical 
measurements such as the gamma-ray and anti- proton data. As a contrast, conventional 
decaying dark matter models have trouble achieving the above goal due to the restrictive 
two-body decay kinematics. 

Using the goldstino as a prime example of the three-body decaying DM, which arises in 
a certain class of supersymmetric theories where SUSY is broken by multiple sectors, we 
found that the goldstino decay with universal leptonic couplings could fit PAMELA and 
Fermi-LAT positron data well, if the goldstino mass is at around 2 TeV and the gravitino, 
which escapes detection and results in missing energy, has a mass around 1 TeV. The mass 
difference is driven by the hardening feature around 300 ~ 500 GeV in the Fermi-LAT 

+ data. A sUghtly smaller (larger) mass difference than 1 TeV can be counterbalanced 
with a shorter (longer) lifetime. 

We further demonstrated that these features of goldstino decays persist in other types of 
three-body decay mechanisms, by studying the injection energy spectra of several different 
types of four-fermi interactions with universal leptonic couplings, as well as scalar dark 
matter decays. We found that these other mechanisms all have softer injection spectra and 
consequently smaller missing particle mass for best fits to the positron data. However, the 
best fit are all similar to the goldstino case. The only exception is when certain operators 
are associated with chiral symmetry breaking and might have couplings proportional to 
the masses of SM fermions. Then the r decay channel would be the dominant one and in 
this situation the resulting soft injection spectrum can not fit e+/ (e~ -|- e"^) and e" -\- at 
the same time, while satisfying the gamma ray constraints. Interestingly, the scalar dark 
matter have the very similar injection spectrum for the mass region of interest and in turn 
the similar best fits. These observations suggest that it would be difficult to distinguish 
different decay mechanisms using data considered in this work. 

One important advantage of three-body decaying dark matter over other conventional 
models is the ability to avoid null searches in cosmic gamma-ray and anti-proton data 
for dark matter. Due to the softer injection energy spectra of the three-body kinematics, 
we showed that the diffuse gamma-ray measurements are compatible with the three-body 
decays, while at the same time maintaining the fits to the PAMELA and Fermi-LAT positron 
data. The allowed hadronic decay widths from the anti-proton data are also larger in the 
three-body decay scenario than in the two-body case, due to the softness of the decay energy 
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spectra. 

In the end, we hope it is clear that there is a strong case for three-body decaying dark mat- 
ter, if the excesses in the positron measurements by PAMELA and Fermi-LAT are beheved 
to be due to dark matter. It would therefore be important to explore ways to definitively 
determine the decay mechanism of the dark matter in other types of measurements, which 
will be a subject for future studies. 
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